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Abstract 

An  example  of  regional  wavefield  analysis  using  three- component  seismic  array  data  is  presented. 
We  examine  the  high-frequency  (0.5-5. 0  Hz)  wavefields  from  three  events  that  occurred  in  or  near 
the  Nevada  Test  Site  (NTS),  two  underground  nuclear  tests  and  one  earthquake,  observed  at  the 
IRIS  Joint  Seismic  Program  three-component  array  located  approximately  400  km  to  the  south  at 
Pinon  Flat  Observatory  (PFO).  The  three-component  array  data  are  first  displayed  as  broadband 
bearing-time  records  (BTR)  which  provide  a  time  history  of  the  propagation  direction  of  coherent 
direct  and  scattered  signals  crossing  the  array;  these  displays  are  used  primarily  as  a  guide  for 
further  analysis.  The  BTRs  reveal  significant  differences  in  the  three  wavefields.  These  differences 
are  somewhat  surprising  in  light  of  the  fact  that  these  data  have  essentially  sampled  the  same 
terrain  between  NTS  and  PFO.  The  BTRs  show  that  the  P  and  S  codas,  and  the  Lg  wavetrains 
for  these  three  events  contain  a  considerable  amount  of  forward  scattered/multipathed  energy.  It 
is  only  with  the  arrival  of  surface  waves  (apparent  velocity  approximately  3  km/s)  that  the  coda 
becomes  more  random  in  nature.  The  fact  that  much  of  this  later  coda  cannot  be  modeled  as  plane 
waves  suggests  that  it  is  not  scattered  from  distant  sources. 

The  Lg  wavetrain  for  one  event  is  examined  in  detail.  A  spectogram,  and  P-SV ,  SH  BTRs 
computed  using  the  three-component  array  data  reveal  that  the  early  portion  of  this  particular 
Lg  wavetrain  is  comprised  of  two  linearly  polarized  SH /Love  waves  and  a  Rayleigh  wave.  The 
most  prominent  arrival  in  this  Lg  wavetrain  is  examined  in  more  detail.  When  plotted  in  three- 
dimensions,  the  tangential  and  elliptical  components  of  particle  motion  are  more  pronounced. 
High-resolution  frequency-wavenumber  analysis  reveals  that  this  section  of  this  particular  Lg  is 
comprised  of  simultaneously  arriving  PSVZ  (principal-component  in  the  P-SV  plane  with  predom¬ 
inantly  vertical  motion)  and  SH /Love  modes  with  different  frequencies,  wavenumbers,  and  slightly 
different  apparent  velocities.  The  nearly  identical  radial  and  Hilbert-transformed  vertical  compo¬ 
nent  waveform  data  provide  further  evidence  that  the  P-SV  energy  is  a  Rayleigh  as  opposed  to  a 
shallow  incidence  angle  SV  wave. 


key  words:  array,  polarization,  seismic  phase  identification,  regional  wave  propagation,  scattering, 
multipathing,  seismic  coda. 
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Objective 


The  objective  of  our  research  is  to:  (i)  improve  event  detection,  phase  identification  and  parameter 
and  waveform  estimation,  (ii)  provide  an  improved  understanding  of  regional  wave  propagation, 
and  (in)  to  aid  in  the  characterization  of  lithospheric  heterogeneities.  These  goals  are  addressed 
by  combining  analysis  of  observed  three-component  array  data  and  numerical  simulations  of  elastic 
wave  propagation  in  realistic  earth  structures. 

•  Preliminary  Research  Results 

For  this  report  we  examine  three  events  that  occurred  in  the  Nevada  Test  Site  area  (NTS);  two 
are  underground  nuclear  tests  and  one  an  earthquake.  These  data  were  collected  as  part  of  the 
Incorporated  Research  Institutions  for  Seismology  (IRIS)  consortium’s  Joint  Seismic  Program  at 
Pinon  Flat,  California  (PFO)  (Al-Shukri  et  al.,  1991).  The  upper  part  of  Figure  1  shows  the  PFO 
array  geometry  and  a  regional  seismicity  map  (from  calnet  and  CALTECH  catalogs,  Ml  >  3)  for 
the  three- month  PFO  experiment.  The  locations  of  NTS  and  PFO  are  shown  on  the  seismicity 
map.  The  lower  part  of  Figure  1  shows  the  broadband  and  0.5-5. 0  Hz  bandpass  filtered  data 
for  the  three  events  examined  for  this  report.  These  waveforms  are  event- wise  normalized  with 
amplitudes  relative  to  the  maximum  in  the  broadband  and  bandpass  filtered  data,  respectively, 
listed  with  each  trace.  Table  1  lists  the  origin  times,  locations,  etc.,  for  the  three  events  examined 
for  this  report. 

•  Processing  Approach 

The  bandpass  filtered  (0.5-5.0  Hz)  three-component  array  data  are  analyzed  using  the  approach 
discussed  in  Wagner  &  Owens  (1995b).  Briefly,  the  magnitude  of  the  multidimensional  enhanced 
minimum- variance  (EMV)  spectrum  is  the  inverse  of  the  minimum  eigenvalue,  A3  (0,ur,c),  of  the 
matrix 

C(0, u,c)  =  E t(0,w,c)  Uy(u/)  Ay  (a;)  uj^w)  E(0,w,c)  (1) 

where  t  denotes  conjugate  transpose,  Uy  (a;)  is  the  unitary  matrix  whose  columns  are  the  [eigen] vectors 
that  span  the  signal  subspace’s  (Owsley,  1978)  orthogonal  complement,  AJy(o?)  is  the  diagonal  ma¬ 
trix  containing  the  associated  eigenvalues,  and 

t  t  t 

E(0,w,e)  =  ez(0,w,c)  eNS(0,w,c)  eEW(0,u;,c)  (2) 

4'  4'  4- 

where  the  et(0,  u,  c)  are  mutually  orthogonal  steering  vectors  for  Z,NS,EW  polarized  plane  waves 
traveling  in  direction  6  with  apparent  velocity  c.  The  eigenvector  associated  with  the  minimum 
eigenvalue,  p3(0,u>,  c),  is  in  general  complex  and  parameterizes  the  component- wise  amplitude  and 
phase  relations.  The  mode  vector  for  the  w  frequency  component  is 

e(0,w,c,p3)  =  E (0,w,c)  p3  (0,  u>,  c)  (3) 

Note  that  Equation  3  merely  provides  high-resolution  estimates  of  the  mode  vectors.  Waveform 
estimates  are  obtained  by  projecting  the  e(0,u>,  c,  p3)  onto  the  signal  subspace  (Gething,  1971) 

w(0, oj, c, p3)  =  A y2(w)  et(0,u7,c,ps)  U5(w)  (4) 
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Figure  1:  Top,  the  PFO  array  geometry  and  regional  seismicity  (Mi  >  3)  that  occurred  during 
the  three-month  PFO  experiment.  Bottom,  the  broadband  and  0.5-5.0  Hz  bandpass  filtered  data 
for  the  three  events  examined  for  this  report.  The  waveforms  are  event-wise  normalized  with 
amplitudes  relative  to  the  maximum  in  the  broadband  and  bandpass  filtered  data,  respectively, 
listed  with  each  trace  See  Table  1  for  origin  times,  locations,  etc.. 
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date 

hr:mn:sec 

lat°  (N) 

Ion0  (E) 

z  (m) 

delta0 

BAZ° 

03/08/91 

21:02:45.40 

37.104 

-116.074 

475 

4.6 

3.49 

5.04 

03/10/91 

21:38:47.51 

37.110 

-115.232 

5,000 

4.0 

3.63 

15.67 

04/04/91 

19:00:00.00 

37.296 

-116.313 

629 

5.4 

3.67 

3.01 

Table  1:  Origin  times,  latitudes,  longitudes,  depths  (meters),  magnitudes,  distances,  and  back 
azimuths  from  PFO  for  the  three  events  examined  in  this  report. 

where  Uj (a>)  is  the  unitary  matrix  whose  columns  are  the  [eigen] vectors  that  span  the  signal  sub¬ 
space,  and  is  the  diagonal  matrix  containing  the  square  roots  of  the  associated  eigenvalues. 

Knowing  the  direction  of  energy  flux,  0,  the  wavetype  can  be  determined  based  on  the  particle 
motion  parameterized  by  p3(0,a>,  c).  The  associated  power  flux  density  is 

V{0,u),c,  p3)  =  et(0,w,c,p3)  Us(u>)  A5(u;)  u|(«)  e(0,u;,c,p3)  .  (5) 

The  following  section  provides  an  outline  of  the  analysis  procedure  using  this  approach. 

•  Data  Analysis 

Our  analysis  typically  begins  by  computing  broadband  bearing-time  records  (BTR)  of  the  three- 
component  array  data  (Wagner  &  Owens,  1995a).  BTRs  provide  a  time  history  of  the  propagation 
direction  of  the  coherent  direct  and  scattered  arrivals  crossing  the  array.  The  BTRs  serve  pri¬ 
marily  as  a  guide  to  further  analysis.  Figure  2  shows  the  broadband  (0.5-5.0  Hz)  EMV  BTRs 
for  the  three  events  shown  in  Figure  1.  Also  calculated  but  not  displayed  in  Figure  2  are  the 
component-wise  amplitude  and  phase  relations  for  the  coherent  modes  (i.e.  the  p3($,  w,  c)).  It 
is  this  information,  combined  with  the  direction-of-arrival  (DOA)  information,  that  allows  us  to 
identify  the  wavetype.  Apparent  velocity  information  is,  unfortunately,  lost  with  this  particular 
broadband  implementation. 

Figure  2  raises  several  questions  of  particular  interest  for  regional  monitoring.  First,  the  in¬ 
creased  interest  in  regional  characterization  is  both  motivated  by  and  based  on  the  assumption 
that  events  that  occur  in  the  same  region  and  sample  the  same  terrain  will  be  similarly  affected 
by  the  structure,  and  that  this  effect  occurs  independent  of  source  magnitude.  One  might  expect, 
therefore,  that  the  wavefields  for  these  three  events  would  exhibit  similar  characteristic.  This  is, 
however,  not  the  case.  The  second  observation  is  that  the  P  and  S  codas,  and  the  Lg  wavetrains 
contain  considerable  amounts  of  forward  scattered/multipathed  energy.  It  is  only  with  the  arrival 
of  the  surface  wave  energy  (apparent  velocity  of  approximately  3  km/s)  that  the  coda  begins  to 
exhibit  more  varied/random  DOAs.  The  fact  that  much  of  this  later  coda  cannot  be  modeled  as 
plane  waves  suggests  that  it  is  not  the  product  of  scattering  from  distant  sources. 

Of  the  three  events  examined  for  this  report,  the  event  that  occurred  on  3/8/91  has  the  best 
developed  Lg  wavetrain.  Because  Lg  is  of  considerable  interest  in  CTBT  monitoring  and  regional 
wave  propagation  in  general,  we  conducted  detailed  analysis  of  the  Lg  wavetrain  for  this  event. 
Figure  3  shows  the  spectogram  (Equation  5  as  a  function  of  time)  for  the  4-3  km/s  apparent  velocity 
window.  The  Z,NS,EW  data  are  plotted  with  the  spectogram  to  help  provide  insight  into  the  way 
each  of  the  components  contributes  to  the  total  field  (note  that  these  data  are  essentially  naturally 
rotated).  Similar  plots  showing  the  relative  magnitudes  of  the  Z,R,T  components  of  p3(9,u,c)  as 
red,  blue,  green  can  be  used  to  help  identify  wavetypes.  Figure  3  shows  the  higher  frequency  Lg 
wavetrain  proceeding  and  superimposed  on  the  lower  frequency  S'- coda  and  surface  wave  energy. 
The  early  Lg  (11-16  secs  in  Figure  3)  is  comprised  of  three  distinct  phases  with  two  simultaneous 
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Figure  2:  Three- component,  broadband  (0.5-5. 0  Hz)  bearing-time  records  for  the  three  events 
listed  in  Table  1. 
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Figure  3:  Spectogram  for  the  4-3  km/s  apparent  velocity  window  for  event  3/8/91. 

arrivals  at  approximately  14.5  secs.  Based  on  p3(0,u;,  c),  these  two  phases  are  identified  as  linearly 
polarized  Love/ SH  and  Rayleigh  waves  with  peak  frequencies  of  approximately  2.5  and  1.8  Hz, 
respectively.  Figure  4  shows  the  BTRs  for  the  total  field,  and  the  field  decomposed  into  P-SV 
and  SH  components  (Wagner  &  Owens,  1993).  These  P-SV,  SH  component  BTRs  reveal  the 
characteristic  simultaneously-arriving  mixed-mode  nature  of  Lg. 

Further  analysis  was  conducted  on  the  most  prominent  arrival  in  the  Lg  wavetrain  for  this 
event  (3/8/91).  The  top  part  of  Figure  5  (part  A)  shows  the  three-dimensional  particle  motion 
for  the  1  sec  window  containing  the  largest  arrival  in  the  Lg  wavetrain.  When  viewed  from  the 
NNW  and  ENE,  the  tangential  and  elliptical  components  of  particle  motion  are  more  pronounced. 
The  middle  part  of  Figure  5  (part  B)  shows  the  EMV  frequency-wavenumber  (u>-k)  spectra  for 
the  PS Vz,  PSVr  and  SH  components  (Wagner  &  Owens,  1993).  The  uo-k  spectra  were  computed 
(using  the  complex  analytic  signal)  for  the  beam  steered  in  the  direction  indicated  by  the  BTRs. 
The  ui-k  analysis  reveals  that  this  data  window  contains  simultaneously  arriving  PSVZ  (principal- 
component  in  P-SV  plane  with  predominantly  vertical  motion)  and  SH  /Love  modes  with  different 
frequencies,  wavenumbers,  and  slightly  different  apparent  velocities.  The  bottom  part  of  Figure  5 
(part  C)  shows  the  vertical-  and  radial-component  waveforms  (left),  and  the  radial-  and  Hilbert- 
transformed  (H)  vertical-component  waveforms  (right).  The  elliptical  particle  motion  (part  A) 
and  the  nearly  identical  H(Z)  and  R  waveforms  concur  with  the  component-wise  phase  relations 
parameterized  by  p3(0,o>,  c),  all  of  which  suggest  that  the  arrival  isolated  on  the  PSVZ  spectrum 
is  a  Rayleigh  as  opposed  to  a  shallow  incidence  angle  SV  wave. 
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Figure  4:  Total  field,  P-SV.  and  SH  bearing-time  records  for  the  0.5-5.0  Hz  passband,  4-3  km/s 
apparent  velocity  window  for  the  event  that  occured  on  3/8/91  (see  Table  1  for  origin  times, 
location,  etc.)- 
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Figure  5:  Top  (A),  three-dimensional  particle  motion,  viewed  from  the  ENE  (left)  and  NNW  (right), 
for  the  1  sec  window  containing  the  most  prominant  arrival  in  the  Lg  wavetrain  for  event  3/8/91. 
Middle  (B),  PSVZ ,  PSVr,  SH  frequency-wavenumber  spectra  for  the  most  prominant  Lg  arrival; 
the  three  spectra  are  plotted  using  the  same  amplitude  sacle.  Bottom  (C),  the  vertical-  and  radial- 
component  waveforms,  and  the  radial-  and  Hilbert-transformed  vertical-component  waveforms  of 
the  early  Lg. 
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•  Conclusions,  Recommendations  &  Future  Plans 

This  report  provides  an  example  of  analysis  of  the  regional  wavefield  using  three- component  seismic 
array  data.  Our  analysis  of  Lg  provides  an  example  of  the  wealth  of  information  provided  by  the 
elastic  wavefield  and  three-component  array  data.  Vertical-component  array  data  provide  no  infor¬ 
mation  about  the  SH / Love  component  of  the  seismic  wavefield,  and  reduce  our  ability  to  detect, 
identify,  and  obtain  truly  representative  waveform  estimate  for  P,  SV  and  Rayleigh  waves.  Elas¬ 
tic  wavefield  decomposition  using  single-site  three- component  data  requires  additional  information 
(i.e.  wavetype,  propagation  direction)  that  may  not  be  available  in  a  realistic  monitoring  situation. 
Three-component  array  data  do  not  restrict  our  analysis  to  the  vertically  polarized  component  of 
the  wavefield,  and  allow  us  to  resolve  the  numerous  problems  and  ambiguities  inherent  in  elastic 
wavefield  decomposition  using  data  form  isolated  three-component  sensors. 

The  complex  problems  we  are  faced  with  in  monitoring  a  CTBT  are  exacerbated  by  not  taking 
full  advantage  of  the  wealth  of  information  provided  by  the  elastic  wavefield.  Elastic  waveform 
data  provide  the  best  and  sometimes  our  only  source  of  information  about  source  processes  and 
earth  structure.  Regional  waveform  data  play  a  vital  role  in  CTBT  monitoring  by  providing  in¬ 
formation  about  smaller  events  that  may  not  be  observed  at  teleseismic  distances.  Broadband 
three-component  array  data  provide  a  multidimensional,  spatially  and  temporally  unaliased  sam¬ 
pling  of  the  regional  wavefield  that  permits  a  nearly  comprehensive  decomposition  of  the  elastic 
wavefield.  These  data  allow  us  to  address  and  answer  fundamental  questions  about  seismic  wave 
propagation  that  simply  cannot  be  answered  conclusively  using  data  from  vertical- component  ar¬ 
rays  or  networks  of  isolated  single-  and/or  three-component  sensors.  The  availability  of  broadband 
three-component  array  data  collected  in  regions  around  the  globe  would  provide  a  unique  and 
invaluable  resource  to  the  entire  community. 

Future  plans  are  to  continue  analysis  of  these  data  to  provide  an  improved  understanding  of 
regional  wave  propagation,  to  improve  event  detection,  phase  identification  and  waveform  esti¬ 
mation,  and  to  aid  in  the  characterization  of  lithospheric  heterogeneities.  Our  analysis  is  limited 
somewhat  by  the  availability  of  three-component  array  data  from  different  regions.  Without  data 
from  other  regions  it  is  not  possible  to  ascertain  which  of  the  characteristics  observed  in  PFO  data 
are  unique  to  this  region  in  particular,  and  which  are  true  for  regional  wave  propagation  in  general. 
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